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A comprehensive analysis of the Symanzik improved anisotropic three-dimensional (7(1) lattice 
gauge theory in the Hamiltonian limit is made. Monte Carlo techniques are used to obtain numerical 
results for the static potential, ratio of the renormalized and bare anisotropies, the string tension, 
lowest glueball masses and the mass ratio. Evidence that rotational symmetry is established more 
accurately for the Symanzik improved anisotropic action is presented. The discretization errors 
in the static potential and the renormalization of the bare anisotropy are found to be only a few 
percent compared to errors of about 20-25% for the unimproved gauge action. Evidence of scaling 
in the string tension, antisymmetric mass gap and the mass ratio is observed in the weak coupling 
region and the behaviour is tested against analytic and numerical results obtained in various other 
Hamiltonian studies of the theory. We find that more accurate determination of the scaling coeffi- 
cients of the string tension and the antisymmetric mass gap has been achieved, and the agreement 
with various other Hamiltonian studies of the theory is excellent. The improved action is found to 
give faster convergence to the continuum limit. Very clear evidence is obtained that in the contin- 
uum limit the glueball ratio Ms /Ma approaches exactly 2, as expected in a theory of free, massive 
bosons. 
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I. INTRODUCTION 

Lattice gauge theory calculations have demonstrated 
important qualitative features of QCD, with increasing 
accuracy. Most lattice gauge theory calculations to date 
have been performed using Monte Carlo techniques in the 
Euclidean formulation. Although the Euclidean lattice 
gauge theory has been a very successful non-perturbative 
technique to compute the properties of elementary par- 
ticles over the years, there are areas where progress has 
been very slow. Examples are QCD at finite temperature, 
glue thermodynamics, heavy quark spectra etc. Some of 
these problems have resisted solution even by the pow- 
erful techniques of Euclidean field theory. This suggests 
that alternative methods should be pursued in parallel 
with Euclidean lattice gauge theory. A viable alternative 
that needs to be explored is the Hamiltonian version of 
QCD. This approach provides a valuable check of the uni- 
versality of the Euclidean results p] and has an appealing 
aspect in reducing lattice gauge theory to a many-body 
problem. As such the formalism is suited for the applica- 
tion of a host of analytic methods imported from quan- 
tum many body theory and condensed matter physics. It 
has been suggested that Hamiltonian lattice gauge theory 
could more readily handle finite density QCD J2J. The 
problems encountered in finite density QCD in the Eu- 
clidean formulation have prompted a return to the strong 
coupling expansions of early Hamiltonian lattice gauge 
theory J3j. Similar ideas have been pursued recently by 
Luo et al., who propose an alternative Hamiltonian 
lattice formulation, 'the Monte Carlo Hamiltonian', and 
have already demonstrated its validity and efficiency for 
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the $ 4 model @. 

Here we attempt to extend the standard Euclidean 
path integral Monte Carlo techniques to Symanzik im- 
proved U(l) gauge theory in three dimensions. Applica- 
tions of this method have been extremely successful Q, Q 
and have given rise to great optimism about the possi- 
bility of obtaining results relevant to continuum physics 
from Monte Carlo simulations of lattice versions of the 
corresponding theory. Such an optimistic view is sup- 
ported by recent work of Sexton et al. 0, Boyd et al. 
[To|. Luo et al. an d Morningstar and Peardon 
who have attempted to derive masses of the low-lying 
hadrons from such calculations and report successful re- 
sults. More recent is the application of PIMC techniques 
to obtain results in the Hamiltonian limit for the U(l) 
model in (2+1) dimensions and SU(3) lattice gau ge 
theory in (3+1) dimensions on anisotropic lattices |14| . 

The use of improved actions [TillTil makes possible ac- 
curate Monte Carlo simulations of QCD on coarse lattices 
with greatly reduced computational effort 0, 0, Il8| . 
In principal, with an improved action it is possible to 
achieve lattice volumes large enough to overcome finite- 
size effects and obtain measurements with good statisti- 
cal errors. Coupled with tadpole improvement the 
pursuit of the Symanzik program has led to significant 
progress in reducing the discretization errors and the 
renormalization of the anisotropy to the level of a few 
percent; makes using anisotropic lattices no more this 
difficult than isotropic ones. At the same time the mer- 
its of using an improved anisotropic lattice have been 
well understood [17], HOj- Anisotropic lattices allow us to 
carry out numerical simulations with a fine temporal res- 
olution while keeping the spatial lattice spacing coarse, 
i.e., a t < a s , where a* and a s are the lattice spacings in 
the temporal and spatial directions, respectively. This is 
especially important for QCD Monte Carlo simulations at 
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finite temperature and heavy particle spectroscopy. But 
more importantly it should make extrapolations to the 
continuum limit more reliable. 

As mentioned above, our aim in this work is to ap- 
ply standard Euclidean path integral Monte Carlo tech- 
niques to extract the Hamiltonian limit for Symanzik im- 
proved C/(l)2+i lattice gauge theory on anisotropic lat- 
tices. The idea is to measure physical quantities on in- 
creasingly anisotropic lattices, and make an extrapola- 
tion to the extreme anisotropic limit. The effect of the 
plaquette improvement is examined by studying the scal- 
ing behaviour and the sensitivity of scaling coefficients of 
the string tension and glueball masses in the Hamiltonian 
limit, £ — * 0. The rest of the paper is organized as fol- 
lows: In Sec. II we briefly review the formulation of the 
Symanzik improved U(l) gauge action in three dimen- 
sions on an anisotropic lattice. The details of the simula- 
tions and the methods used to extract the observables are 
described in Sec. III. Here we discuss our techniques for 
calculating the static quark potential, renormalization of 
anisotropy, string tension and glueball masses from Wil- 
son loop operators. We present and discuss our results 
in Sec. IV. Scaling of the string tension, antisymmetric 
mass gap and the mass ratio in the weak-coupling re- 
gion, are tested against theoretical predictions and com- 
pared with the estimates obtained by other studies in the 
Hamiltonian limit. Our conclusions are given in Sec. V, 
along with an outline of future work. 



II. IMPROVED ANISOTROPIC 
DISCRETIZATION OF U(l) 2 +i 

The Symanzik improved U (1)2+1 gauge action on an 
anisotropic lattice is identical in form to the SU(3) case 
and is given 1 by 01 
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where P„„ and are the lxl Wilson loop and 2x1 
rectangular loop in the /i x v plane respectively. At the 
tree-level the coefficients are chosen so that the action 
has no 0(a 2 ) discretization corrections. The spatial and 
temporal square and rectangular loops are given by 



Pij(x) = l-U i (x)U j (x + i)Uj(x + j)W(x) 



(2) 



Rij(x) = l-Ui(x)Ui(x + i)Uj(x + 2i)Uj(x + i + j) 



1 The notation used here differs slightiy from that used in Ref. 
|2l| , where the prefactors were absorbed into the definitions of /3 
and £. We follow the notation introduced in Il7l . 
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1 - Ui(x)U t (x + i)Uj(x + t)ul(x) 
1 - Ui(x)Ui(x + i)U t (x + 2i)JjJ(x + i + t) 
xUj(x + t)U}(x)], (5) 



where x labels the sites of the lattice, i, j are the spatial 
indices and Ui(x) is the link variable from site x to x + i. 
The rectangular loop that extends two steps in the time 
direction has not been included in the above action. This 
has the advantage of eliminating negative residue high 
energy poles in the gluon propagator 17] but at the same 
time leaves errors of the order a 2 in the action. The 
0(a 2 ) errors, however, are negligible provided at is small 
compared to a. The bare anisotropy parameter is £ and 
is equal to the aspect ratio of the temporal and spatial 
lattice spacings at the tree level. At higher orders in 
the perturbative expansion, the bare anisotropy in the 
simulated action is not the same as the measured value, 
£phys, due to quantum corrections. The couplings (3 S and 
(3t are defined by 
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The two different couplings in Eq. are necessary in 
order to ensure that in the continuum limit, physical ob- 
servables become independent of the kind of lattice reg- 
ularization chosen. In the case of an asymmetric lattice, 
this implies that physical quantities have to be indepen- 
dent of the anisotropy factor £. To achieve this one needs 
to introduce different couplings for spatial and temporal 
directions, which depend on £. The ^-dependence of the 
couplings g 2 and g\ is due to quantum corrections and 
leads to an energy sum rule for the quark-antiquark po- 
tential, and the glueball mass, which differs in an impor- 
tant way from that which one would expect naively. 

In the weak coupling limit of SU(N) lattice gauge 
theory, the relation between the scales of the cou- 
plings in Euclidean and Hamiltonian formulations has 
been determined analytically from the effective actions 
Ullll HI, using the background field method on 
the lattice. For small values of gE, the couplings g s and 
gt can be expanded as 
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(7) 
(8) 



where gE is the Euclidean coupling. For £ = 1, one re- 
covers the usual Euclidean lattice gauge theory, where 
g s = g t = gE- In the limit £ —> 0, Eqs. (7) and 
(8) reduce to their Hamiltonian values and one obtains 
the relation between the Euclidean coupling gE and its 
Hamiltonian counterpart gu- Similar calculations have 
been performed to determine the anisotropic coefficients 
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c s and ct at arbitrary anisotropy for a class of improved 
actions [2(], H^, HU. These coefficients become an 
imp ortant tool in the analysis of glue thermodynamics 
|29| , the quark-gluon plasma [3(]j El and for the deter- 
mination of spectral functions at finite temperature |32| . 
Similar calculations have not yet been done for the U(l) 
theory, however, as far as we are aware. 

Tadpole improvement [Tflj is introduced by renormal- 
izing the link variables: Ui(x) — > Ui(x)/u s , and Ut(x) — > 
/7 t (x) / Ut , where the mean fields u s and u t can be defined 
by using the measured values of the average plaquettes 
in a simulation. In the plaqucttc mean-link formulation, 
the mean fields are determined self-consistently and are 
defined by 

U, = (Py) 1/4 , U 2 t u\=P lt . (9) 

For at << a, the mean temporal link ut is expected to 
be very close to unity. For simplicity we use a convenient 
and gauge invariant definition for u s in terms of the mean 
spatial plaquette given by u s = (ReP^) 1 / 4 , and compute 
ut from the temporal plaquette Pa in Eq. ©. 

To obtain the Hamiltonian estimates from anisotropic 
lattices, a naive extrapolation procedure is followed. In 
this procedure we assume classical values of the cou- 
plings, i.e., /3 = p s = fa in Eq. Q and extrapolate 
the physical quantities to the extreme anisotropic limit, 
£ — > at constant f3. Such a procedure is not strictly cor- 
rect, however, at the quantum level because /3 S ^ fit 7^ /? 
due to renormalization 2 . 

As an example of the application of PIMC to the 
Symanzik improved action, we consider the case of com- 
pact 17(1) gauge theory in three dimensions. The rele- 
vance of the model to QCD at finite temperature [34| has 
made it a standard proving ground for Hamiltonian lat- 
tice numerical methods. The model has two essential fea- 
tures in common with QCD, confinement [351I361 13H l3S| 
and chiral symmetry breaking |39|. Other common fea- 
tures are the existence of a mass gap and a confinement- 
deconfinement phase transition at some non-zero tem- 
perature. These similarities suggest that a comparison 
of the respective mass spectra should be informative. In 
the continuum limit of theory, the mass gap M is found 
to behave as |36| 

AfV « const./3exp[-27r 2 /3t;(0)] = M 2 D (10) 

where v(0) « 0.2527 is the lattice Coulomb Green's func- 
tion at zero separation. In Hamiltonian theory in which 
the space dimensions are discretized, vh(0) ~ 0.3214 is 
the analogous Green's function [33. It has been shown 
analytically that a linear confinement persists for all non- 
vanishing couplings, no matter how weak [3(|[33]. The 
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and the gauge coupling in spatial and temporal directions for 
the improved Abelian lattice gauge theory are currently under 
way and will be reported elsewhere l33l . 



string tension as a function of coupling also scales expo- 
nentially and obeys a lower bound 

K > const.Mn/T 1 . (11) 

An interesting feature to explore in this context is 
whether the coupling to matter fields will change the per- 
manent confinement status in (2+1) dimensions ppf. 

III. METHOD 
A. Simulation details 

To extract estimates of the string tension and the glue- 
ball masses using the Symanzik improved action in Eq. 
(fT|). a set of simulations are performed on lattices of size 
TV 2 x N t , (N s = 16, and N t = 16 - 64) where N s and 
Nt are the number of lattice sites in the spatial and tem- 
poral directions respectively. The lattice size in the time 
direction is adjusted according to the anisotropy used in 
order to keep the physical length in the spatial and tem- 
poral directions equal. To analyse the behaviour in the 
strong and weak coupling regions, gauge configurations 
are generated using the Metropolis algorithm, for a range 
of couplings P = 1 — 2.5. The details of the algorithm 
are discussed elsewhere |l3j . Starting from an arbitrary 
initial gauge configuration, 50,000 sweeps are performed 
for the equilibration of the configurations and the self- 
consistent determination of the mean-field parameters. 
A Fourier acceleration procedure 0, 0] ^ s use d to over- 
come the stiffness against variations in the temporal pla- 
quettes for high anisotropics. About 50% of the ordinary 
Metropolis updates are replaced with Fourier updates for 
£ < 0.4 and about 100,000 further sweeps are performed 
to allow the system to equilibrate. 

After thermalization, configurations are stored every 
300 sweeps; 1200 stored gauge configurations are used in 
the measurement of the static quark potential and string 
tension and 1500 configurations for the glueball masses at 
each [3. Measurements made on the stored configurations 
are binned into five blocks with each block containing an 
average of 250 measurements. The mean and standard 
deviation of the final observables are estimated simply 
by averaging over the block averages. The simulation 
parameters used for each configuration set are shown in 
Tabled 

B. The inter-quark potential and string tension 

The static quark potential V(r) is extracted from the 
expectation values of the time-like Wilson loops W(r, t), 
which are expected to behave as 

iy(r,t)«^Z 4 (r)exp[-^(r)], (12) 

i 

where the summation runs over the excited state contri- 
bution to the expectation value, and i = 1 corresponds 



TABLE I: Simulation parameters at various f3 and £ vales. 



Volume 


P 




u t 


U 3 


/ \4 

M 


< P > at it s 


16 2 


x 32 


1.0 


0.5 


0.9991(3) 0.9218(2) 


0.7220 


0.7221(4) 


16 2 


x 32 


1.45 


0.5 


0.9991(2) 0.9468(2) 


0.8038 


0.8039(3) 


16 2 


x 32 


1.75 


0.5 


0.9997(3) 0.9551(3) 


0.8324 


0.8325(6) 


16 2 


x 32 


2.0 


0.5 


1. 


0.9589(2) 


0.8458 


0.8459(4) 


16 2 


x 32 


2.5 


0.5 


1. 


0.9614(4) 


0.8543 


0.8544(3) 


16 2 


x 40 


1.0 


0.4 


1. 


0.9182(3) 


0.7110 


0.7111(4) 


16 2 


x 40 


1.45 


0.4 


1. 


0.9432(3) 


0.7917 


0.7918(7) 


16 2 


x 40 


1.75 


0.4 


1. 


0.9497(3) 


0.8138 


0.8139(5) 


16 2 


x 40 


2.0 


0.4 


1. 


0.9524(5) 


0.8228 


0.8229(5) 


16 2 


x 40 


2.5 


0.4 


1. 


0.9561(2) 


0.8358 


0.8359(2) 


16 2 


x 48 


1.0 


0.333 


1. 


0.9152(3) 


0.7016 


0.7017(5) 


16 2 


x 48 


1.45 


0.333 


1. 


0.9356(3) 


0.7663 


0.7664(4) 


16 2 


x 48 


1.75 


0.333 


1. 


0.9401(2) 


0.7815 


0.7816(6) 


16 2 


x 48 


2.0 


0.333 


1. 


0.9408(4) 


0.7836 


0.7837(4) 


16 2 


x 48 


2.5 


0.333 


1. 


0.9454(4) 


0.7991 


0.7992(3) 


16 2 


x 64 


1.0 


0.25 


1. 


0.9112(6) 


0.6895 


0.6892(3) 


16 2 


x 64 


1.45 


0.25 


1. 


0.9287(3) 


0.7440 


0.7439(4) 


16 2 


x 64 


1.75 


0.25 


1. 


0.9357(5) 


0.7665 


0.7663(5) 


16 2 


x 64 


2.0 


0.25 


1. 


0.9381(4) 


0.7745 


0.7738(6) 


16 2 


x 64 


2.5 


0.25 


1. 


0.9401(4) 


0.7812 


0.7811(3) 



the logarithmic ratio of successive Wilson loops 

"W(r,i + 1) 



Vt{r) 



-In 



W(r,t) 



(14) 



which is expected to be independent of t for t > 0. A plot 
of the effective potential is shown in Fig. [I] for (3 = 1.70 
and £ = 0.4 for various separations r. The dashed lines 
indicate the plateau values at various separations. As 
a result of heavy smearing, a good plateau behaviour is 
seen at small t values for r = 1 through 7. For r > 6, we 
see that the signal is dominated by noise for t > 6. We 
fix the fitting range to be, in most cases, t = 2 to 6. 
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to the lowest energy contribution. To obtain the optimal 
overlap of Wilson loop (and glueball) operators with the 
lowest state, it is necessary to suppress the contamina- 
tion from excited states. This is done by using a simple 
APE smearing method 0, 0, El which is implemented 
by the iterative replacement of the original spatial link 
variables by a smeared link. Following the single-link 
smearing procedure, every space- like link variable Ui(x) 
on the lattice is replaced by 



Ui(x) -> P u{1) 



where [7_j = U'(a 



aUi{x) U j( x ) U i( x ± 3) U] ( X ± *) 

(13) ■ 

i) and i and j are purely spatial 



indices. Pu(i) denotes the projection onto U(l) and a 
is the smearing parameter. Operators constructed out 
of smeared links dramatically reduced the mixing with 
high frequency modes, thus removing the excited-state 
contamination in the correlation functions. The smearing 
fraction is fixed to a — 0.7 and ten iterations of the 
smearing process are used. To reduce the variance, we 
use the technique of thermal averaging Il3l |45| . which 
amounts to replacing the time-like link variables Ut by 
their local averages. This technique was applied to all 
temporal links except those adjacent to the spatial legs 
of loops, which are not independent. The technique has 
a dramatic effect in reducing the statistical noise. 

The values of the effective potential are measured from 



FIG. 1: Effective potential as a function of Euclidean time 
t. From bottom up the horizantal lines correspond to r = 1 
through 7 and indicate the plateau values in the range 2 < 
t < 6. 

The string tension is then extracted from the Wilson 
loops by establishing the linear behaviour for the static 
quark potential at large separation r. We have chosen to 
fit our results for V(r) to the form [37j 



V(v) 



Kr + Wn(r), 



(15) 



where the linear term dominates the behaviour at large 
separations and a logarithmic Coulomb term dominates 
at small separations. 



C. Renormalization of Anisotropy 

Since the anisotropy ratio, r), is important in QCD 
simulations on the anisotropic lattices, we study its be- 
haviour by numerical simulations. Measurements of 
renormalization of anisotropy |46j have been made by 
comparing the static quark potential extracted from cor- 
relations along the different spatial and temporal direc- 
tions. On an anisotropic lattice there are two different 
potentials Vt(r) and V s (r) extracted from two different 
types of loops: time-like W x t and space-like W xy Wilson 
loops. The two potentials differ by a factor of (, p hys and 
by an additive constant, since the self-energy corrections 
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to the static quark potential are different if the quark and 
anti-quark propagate along the temporal or a spatial di- 
rection. The natural way to proceed is to build ratios of 
the Wilson loops. 

W xt (x,t+1) 



Rt(x,t) = 

R s {x,y) = 

Asymptotically for large r 
approach 

,-tV. , 



W xt (x,t) 
W xy (x,y + 1) 
W xy {x,y) 1 ° j 

and y, the ratios Rt and R s 



Rt{x,t) 
R s (x,y) 



Z XT G 



,-yV s 



(excited state contr.) 
(excited state contr.) 



(17) 



The physical anisotropy is then determined from the ratio 
of the potentials Vt(r) and V a (r) estimated from Rt and 
R s respectively. The lattice potentials defined by Eq. 
(|17|l contain contributions from the self-energy terms. 
The potential is simply parameterized as 

V.(t,r) = V?®+Vf(£,r), (18) 

where V/ is the lattice potential free of self-energy con- 
tributions. The time- like potential V* is treated similarly. 
To eliminate the effect of the self-energy term V° in the 
potentials, we define a subtracted potential 



V f "*(t,t,to) = V/(t,t)-V/(t,to) 



(19) 

where the subtraction points r$ and to are chosen to 
satisfy to = £ro and the matching of the potential 
Vt (*o = £ r o) = V/(ro) should be satisfied there. The 
subtraction radius ro should be chosen to be as small as 
possible so that fluctuations of the potential does not in- 
crease in which case simulations with high statistics on a 
larger lattice are required. The renormalized anisotropy 
is determined from the ratio 



£phys 



ysub 



(20) 



Vf*&r, r ) 

and the measured anisotropy parameter, 77, is then given 
by 

£,phys 



V 



(21) 



An alternative approach to make the comparison is to fit 
the measured potentials with the forms |47| 

a s V s (x) = a + <ra 2 s x + c\nx 

&sVt(t) = a + cra s a t t + cln — — (22) 

a t t 

The renormalized anisotropy £, p hy S is then determined 
from the ratio of the coefficients of the linear terms in 
the two cases. It is advantageous to use the potential 
at smaller r, where the statistical errors are smaller, and 
determine t; p hys from the ratio of the coefficients of the 
Coulomb terms; however, such an estimate depends on 
short distance effects and is more sensitive to possible 
discretization errors of 0(a 4 /r 4 ) |4t| . 



D. Glueball masses 

The numerical analysis of the mass of a glueball hav- 
ing a given J PC proceeds through a study of the time- 
like correlations between space- like Wilson loop operators 

*<(*) 



where 



at) = mmm, 



*i(t) = $<(*) - (0 1 *<(t) 1 0) 



(23) 



(24) 



is a gauge invariant vacuum subtracted operator capable 
of creating a glueball out of the vacuum. It is necessary to 
to subtract the vacuum contribution for the scalar glue- 
ball with J PC = ++ , because in the large Euclidean 
time limit, C(t) becomes dominated by the lowest en- 
ergy state carrying the quantum numbers of $ and these 
quantum numbers may coincide with that of the vacuum. 
The vacuum contribution is averaged over the whole en- 
semble before subtracting from the correlator. The glue- 
ball mass of interest is then extracted by studying the 
exponential decay of the correlation function for large 
Euclidean times, which is expected to behave as 



C(t) = Ci [exp(— rriit) + exp(— m^(T — t))] 
+ (excited state contributions) 



(25) 



where m, is the mass of the lowest-lying glueball which 
can be created by <&i(t), and T = N t at is the extent of 
the periodic lattice in the time direction. Here, only the 
lowest "symmetric" (PC = ++) and "antisymmetric" 

(PC = ) glueball states are studied. The measured 

values of C (t) are expected to fall on a simple exponential 
curve assuming that the lattice is fine enough for the 
glueball mass to exhibit scaling behaviour according to 
the theoretical predictions. 

On a finite lattice with lattice spacing a, the operator 
<&i(t) has a small overlap with the glueball ground state 
and the mass extracted from C(t) may be too large ow- 
ing to the excited-state contamination. The overlap gets 
worse as the lattice spacing is reduced and nears the con- 
tinuum limit, a — > 0. This is obvious since the physical 
extension of the glueball remains fixed, while the operator 
$, constructed from small loops, probes an ever smaller 
region of the glueball wave function as the lattice spacing 
is decreased. Hence it becomes important to use an im- 
proved glueball operator so as to have approximately the 
same size as the physical size of the glueball. For such 
an operator, the overlap with the glueball of interest is 
strong at small lattice spacing and signal-to-noise ratio 
is also optimal |48| . 

Following the variational technique of Morningstar and 
Peardon and the smearing procedure of Teper [49) . 
an optimized operator was found as a linear combination 
of the basic operators <j> 



(26) 
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where the index a runs over the rectangular Wilson loops 
with dimensions l x = [n — l,n + 1], l y — [n — l,n + 1] 
and smearing n s = [m — 1, m + 1]. The correlation func- 
tion C(i) is then computed from the optimized glueball 
operator <&i(t) 



a(t) = J2io I *i(t + *o)$<(*o) I 0) 



(27) 



Fig. |2]shows the effective mass plot for the scalar glue- 
ball for the measurements at /3 = 2.0 and £ = 0.50. The 
signal is seen out to the time slice 6 and reaches a plateau 
region for 1 < t < 5. The data are noisy for t > 5. 




FIG. 2: Effective mass plot for the scalar glueball at (3 — 2.0 
and £ = 0.50. 

We choose to fit the correlation function C(t) with the 
simple form 

C, (t) = a cosh mi [T/2 - 1] (28) 
to determine the glueball mass estimates. 

IV. RESULTS AND DISCUSSION 



Y(r> 




FIG. 3: Static potential, V(r), as a function of separation r 
at j3 = 1.55 and £ = 0.4. The solid line is a fit to the form 
V(r) — a + br + cln(r) in the range 3 < r < 8. 



V(r>3 




FIG. 4: A close up of the static quark potential, V(r), at 
small r, using the Symanzik improved action. The solid line 
is a fit of Eq. 11151 to on-axis points r = 4 to 8. This plot in- 
volves measurements at (3 = 1.55 for £ = 0.4 with 10 smearing 
sweeps at smearing parameter a = 0.7 



A. Static quark potential and rotational symmetry 

A plot of the static quark potential V(r) as a function 
of r at /3 = 1.55 and £ = 0.4 is shown in Fig. The 
data in this plot were obtained by looking for a plateau 
in the effective potential. 

Because we are concerned to make long distance be- 
haviour consistent in both fine and coarse directions, it 
is advantageous to use Wilson loops of the largest pos- 
sible spatial extent. However, in practice, the statistical 
errors in large Wilson loops grow exponentially with sep- 
aration. We use Wilson loops of size 8x8 and fit the data 
by H15fl . We fit in the range 3 < r < 8 so that we are 
not sensitive to the Coulomb term and the discretization 
errors associated with it. We see that the data are fitted 
very well giving the string tension, K(= era 2 ) = 0.044(3). 



One of the main features of the improved discretization 
is the improved rotational invariance [211 ] . Discretization 
errors in the gluon action affect the extent to which con- 
tinuum symmetries, such as rotational symmetry, are re- 
stored. To explore the extent to which rotational invari- 
ance is improved, we measure the potential at off-axis as 
well as on-axis separations. Thus for improved rotational 
invariance, the static potential, for example, at r — (4, 3) 
should agree exactly with that at r = (5, 0). 

For a fixed number of configurations and constant 
physical volume, we show the results from the Symanzik 
improved action and standard Wilson action in Figs. 01 
and[5]respectively. We see that the off axis points for the 
improved lattice are excellently fitted by the rotationally 
invariant fitting curve (0 through r = 3 to 8. The data 
from the Wilson action lie rather less close to the line 
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of best fit. However, at large separations, the standard 
Wilson action does just as well as the improved action in 
extracting the static quark potential. 



V(r) 




FIG. 5: A close up of the static quark potential, V(r), at 
small r. The solid line is a fit of Eq. 11511 to on-axis points r — 
4 to 8. This plot involves measurements from the standard 
Wilson action [lj| at /3 — 1.4142 for £ = 0.4 with 10 smearing 
sweeps at smearing parameter a — 0.7. 

As a quantitive measurement of the improvement, the 
potential measured in the simulation from non-planar 
Wilson loops is compared with an interpolation to the 
on-axis data lial: 



AV(r) 



V si m(r)-V &t (r) 



(29) 



Results for r = (1,1) are given in Table ITU With the 
mean- field inspired Symanzik improvement, the differ- 
ence is only a few percent compared to a difference of 
about 10-20% for the Wilson action |T^ . 



B. Numerical determination of renormalized 
anisotropy 

We choose the points ro = 2 and \/2 to compute the 
subtraction potentials of Eq. (|19JI an d use them to ob- 
tain the ratio in Eq. I|2(J|I . The subtracted spatial and 
temporal potentials at the subtraction point ro = V2 are 
shown in Fig. [5] The anisotropics measured at these sub- 
traction points for the improved and unimproved actions 
at different (3 values are compared in Table ITTI The two 
determinations of the anisotropies are in excellent agree- 
ment. These results show that the input anisotropy is 
normalized by less than a few percent for the improved 
action. This is in contrast with the standard Wilson ac- 
tion |l3j |. where the measured anisotropy is found to be 
about 20-30% lower than the input anisotropy £. 

Fig. shows the potentials computed from spatial 
and temporal Wilson loops without subtracting the self 
energy terms. We find that the difference between the 
estimates computed from subtracted and unsubtracted 



TABLE II: Comparison between the measured anisotropy 
£,phys and the input anisotropy £ for the Symanzik improved 
and the standard Wilson actions |T^I . The measured differ- 
ence in the off-axis potential at r = (1, 1) are also shown. 



Action 


P 




£phys 

r () = 2 ro = y/2 


AV(V2) 


Improved 


1.35 


0.50 


U.4yO^D ) 




0.03(1) 


action 




0.40 




U.oyU^D ) 


0.03(2) 






0.333 


U.ozo^o ) 


U.ozU^o ) 


0.04(4) 




1.45 


0.444 




n A A 1 (A \ 


0.04(2) 






0.333 


U.oZo{o ) 




0.04(3) 






0.25 


0.246(4) 


0.240(7) 


0.05(4) 




1.55 


0.40 


0.402(2) 


0.398(6) 


0.04(2) 






0.25 


0.241(5) 


0.239(7) 


0.06(4) 




1.65 


0.333 


0.334(5) 


0.332(5) 


0.02(1) 




1.75 


0.25 


0.252(4) 


0.249(6) 


0.03(1) 




2.0 


0.333 


0.335(6) 


0.334(2) 


0.03(1) 






0.25 


0.242(4) 


0.246(1) 


0.06(2) 


Wilson 


1.35 


0.444 


0.418(5) 


0.419(7) 


0.12(1) 


action 


1.55 


0.40 


0.379(3) 


0.374(6) 


0.10(1) 




1.70 


0.333 


0.286(7) 


0.289(5) 


0.08(1) 




2.0 


0.333 


0.281(6) 


0.292(2) 


0.13(2) 
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FIG. 6: Subtracted spatial and temporal potentials against 
the separation r. This plot involves measurements at f3 = 1.75 
and f = 0.333 for the subtraction point ro = %/2. The tempo- 
ral potential Vt ub has been rescaled by the input anisotropy. 



potentials is less than 1% for the Symanzik improved ac- 
tion. We conclude that a few percent renormalization 
in the anisotropy is sufficiently small that it is unlikely 
to represent the dominant effect in the final estimates 
and it is safe to use the bare anisotropy for the tadpole 
improved Symanzik action. 
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V(r) 




FIG. 7: Unsubtracted spatial and temporal potentials against 
the separation r at /3 = 1.75 and £ = 0.333. The temporal 
potential Vt(r) has been rescaled by the input anisotropy. 



C. String tension 

To obtain estimates of the string tension in the Hamil- 
tonian limit, an extrapolation is performed by a simple 
quadratic fit in powers of £ 2 for each (3 value. The simu- 
lations run over a range of anisotropies, £ = 1 — 0.25, thus 
enabling reliable extrapolation to the Hamiltonian limit. 
The errors for the extrapolation may be obtained by the 
"linear, quadratic, cubic" extrapolant method [5(J • Fig. 
[5] shows our estimates of the string tension as a function 
of the anisotropy £ 2 for various fixed (3 values. Except 
at (3 = 1.35, a fairly smooth variation of string tension 
with £ 2 for various couplings is seen. The curvature in 
the extrapolation at f3 — 1.35 suggests that our estimate 
may be somewhat too high there. 




FIG. 8: Extrapolation of the string tension to the Hamilto- 
nian limit £ — *• 0, for various f3. Solid lines show the quadratic 
fits in £ 2 to the data. 

Our extrapolated results for the string tension, K = 
a 2 cr, together with the earlier Hamiltonian estimates ob- 
tained from the t-expansion [Hl| , Green's Function Monte 



Carlo simulation s 1531 and the Exact Linked Cluster Ex- 
pansion (ELCE) [53| are plotted as a function of inverse 
coupling in Fig. El We see that the string tension displays 
an exponential behaviour at weak coupling in accordance 
with the theoretical prediction. It has rigorously been 
shown that the string tension in ?7(l)( 2 +i) undergoes a 
'roughening transition 1 [H3 . l55j | at some intermediate cou- 
pling estimated to be near f3 « 0.8. Beyond the transi- 
tion point, the different estimates of the string tension 
are expected to agree. The on-axis strong coupling se- 
ries approximants fail to converge beyond (3 = 0.8, which 
prevents the analytic continuation of the series expan- 
sion beyond the roughening transition. The i-expansion 
results, however, do not suffer from this difficulty. A com- 
parison with the GFMC [52| and an exact linked-cluster 
expansion 53] and results shows that our estimates are in 
good agreement with earlier estimates. The i-expansion 
estimates [Hi] are a little high, but still reasonably accu- 
rate. Of course, we do not expect that the results for the 
improved action should match exactly at finite coupling 
with other estimates which were computed for the unim- 
proved action. We believe our PIMC estimates are more 
reliable and accurate and are also clearly consistent with 
the behaviour predicted by Polyakov |3y| and Gopfert 
and Mack 



Improved action ^A-i 

Wilson action ^ 1- 

t expansion " •- 

ELCE - a 

GFMC ^-O^ 



1 1.2 1.4 1.6 



2.2 2.4 



FIG. 9: The string tension as a function of inverse coupling. 
Our present estimates are shown as open triangles. Earlier 
results from standard Wilson action |l3j|. i-expansion |Hl|. 
Greens Function Monte Carlo simulations |5^| and an exact 
linked cluster expansion (53^1 are shown as solid triangles, solid 
circles, open circles and open squares respectively. 

Fig. 1101 shows the scaling behaviour of the string ten- 
sion together with results obtained using the standard 
Wilson action as a function of (3. The dashed- 

dot line is the strong-coupling expansion to order (3 A 
|5t| and the solid line represents a fit to the weak- 
coupling asymptotic form ©. It can be seen that our 
present estimates appear to match nicely onto the strong 
and weak coupling expansions in their respective limits. 
The 4th order strong coupling series expansion, obtained 
from integrated differential approximants, diverges be- 
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yond (3 = 1.30. In the weak-coupling region, the string 
tension is consistent with the predicted scaling behaviour 
[8?j | . An unconstrained fit of the form © represents the 
data rather well in the interval 1.45 < (3 < 1.95. The 
fit to the data gives a scaling slope of 3.04 ± 0.13 and 
an intercept of 1.42 ± 0.21. The intercept of the scaling 
curve is roughly two times larger in magnitude than the 
theory predicts, compared to our previous results with 
the standard Wilson action ]t| which were higher than 
theory by a factor of 5 - 6. Also in contrast with the Wil- 
son action, a significant reduction in the errors is clearly 
apparent with the tadpole-improved Symanzik action. 

In summary, it appears that the overall exponential 
scaling behaviour is the same for both actions, but the 
constant coefficient is lower for the Symanzik action by 
a factor of 2 to 3, and closer to the theoretical weak- 
coupling estimate. It seems highly plausible that a dif- 
ferent action should give a different renormalization for 
the constant coefficient, although no analytic calculation 
of this effect has been done. 







Improved action ' A ' 
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1 1.2 1.4 1.6 1.8 2 2.2 2.4 2.6 

P 

FIG. 10: The logarithm of the string tension as a function 
of inverse coupling. Our present estimates are shown as open 
triangles. The solid line is the result of fitting to form © for 
1.45 < (3 < 1.95. The dash-dot line is the f3 4 order strong cou- 
pling expansion [5ri||. Our previous estimates flU are shown 
as solid triangles. The dashed line represents the finite size 
behaviour 15711 . 



D. Antisymmetric mass gap 

The weak-coupling behaviour of the mass gap is not ex- 
actly known. The rigorous analysis of Gopfert and Mack 
|37| showed that in the continuum limit [/(l) 2 +i reduces 
to a massive scalar free field theory, with a mass gap M 
which is expected to decrease exponentially as the lat- 
tice spacing goes to zero. They showed that the lattice 
photon mass in the Villain action on a 3-dimcnsional Eu- 
clidean lattice is given by Eq. ijfJJl. It is often claimed in 
the literature that the Villain action is a high-/? approx- 
imation of the Wilson action so that Eq. © should also 



hold in the weak-coupling limit of the Wilson model. 

The extrapolation of the glueball masses to the Hamil- 
tonian limit is shown in Fig. ^] The extrapolation is 
performed by using a simple quadratic fit in powers of 
£ 2 . Again we see a smooth dependence on £ 2 , for all (3 
values analysed here. 
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FIG. 11: Extrapolation of the antisymmetric glueball mass to 
the Hamiltonian limit £ — » 0. Dashed lines show the quadratic 
fits in £ 2 to the data. 

Our extrapolated results for the mass gap together 
with earlier Hamiltonian estimates obtained from the 
Wilson action [Tj| are plotted as a function of (3 in Fig. 
1121 Unlike the string tension, the strong-coupling ex- 
pansion of the mass gaps is believed to be analytic near 
the roughening point. Comparison with earlier Hamilto- 
nian estimates shows that our present data follow quite 
closely the strong-coupling expansion estimates, obtained 
by the method of integrated differential approximants, in 
the strong and weak-coupling regions. The i-expansion 
estimates [Hi], obtained from D-Pade approximants, are 
found to be substantially less accurate than those from 
the SC-expansion. Beyond (3 — 2, the approximants for 
the t-expansion do not converge well so that no reliable 
estimates can be obtained beyond that coupling value. 

TABLE III: Comparison of the Hamiltonian results obtained 
for the symmetric and antisymmetric scalar glueball masses 
(Ms, Ma) for the Symanzik improved and Wilson actions. 



f3 Improved action Wilson action 





M s 


M A 


M s 


M A 


1.0 


1.68(1) 


1.510(5) 


1.9(1) 


1.45(7) 


1.35 


1.207(9) 


0.811(5) 


1.4(1) 


0.9(1) 


1.55 


0.75(1) 


0.444(5) 


0.8(2) 


0.4(1) 


1.70 


0.69(2) 


0.361(6) 


0.5(3) 


0.24(9) 


1.90 


0.352(6) 


0.177(3) 


0.3(1) 


0.17(9) 


2.0 


0.272(4) 


0.136(2) 


0.2(1) 


0.10(6) 



The asymptotic scaling behaviour of the antisymmetric 
mass gap is shown in Fig. The solid line is the result 
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FIG. 12: Glueball masses as functions of /3. Our present 
results for the symmetric and antisymmetric glueball masses 
(Ms, Ma) are shown as open circles and triangles respectively. 
Earlier results obtained from the Wilson action fl3| are shown 
as solid circles and triangles respectively. 



FIG. 13: The logarithm of the antisymmetric mass gap as 
a function of /3. Open triangles show our present estimates. 
The solid line is the fit to the data for 1 < j3 < 2. Previ- 
ous estimates from the Wilson action fl^ are shown as solid 
triangles. 



of fitting for 1.4 < (3 < 2.25 to the form (0 to find the 
scaling slope and the intercept of the scaling curve. Our 
results for these coefficients are shown and compared with 
previous studies in Table llVl These results are obtained 
by fitting to the form M 2 = /3exp[— /o/3 + /i] in the weak- 
coupling region. It can be seen that the agreement with 
the earlier results is remarkable. We find the constant co- 
efficient (intercept of the scaling curve) is approximately 
1.5 times larger in magnitude than the theory predicts. 
This is an improvement over our previous estimate using 
the Wilson action 13], where the constant coefficients 
were estimated a factor 5-6 times larger. The scaling 
slope is a little less than the theoretical prediction but in 
agreement with the estimates obtained in other numerical 
and analytic calculations (Ref. Table llV|) . Several stud- 
ies have provided evidence that the antisymmetric mass 
gap in the Wilson model of t/(l)2+i Hamiltonian lattice 
gauge theory does not fall in the weak-coupling limit in 
the same manner as the periodic Gaussian model 3 . This 
may be a signal of non- universality for the Abclian U(l) 
theory. 



E. Mass gap ratio 

The quantity of interest here is the dimensionless ratio 
Rm between the symmetric and antisymmetric mass gap 
in the large (3 limit. In this limit the ratio Rm = Ms /Ma 
is expected to tend smoothly to its continuum value. In 



3 Based on the fact that periodic Gaussian mod els are special forms 
of Wannier-fuction expansions, Suranyi l58l has argued that a 
natural series of models, beginning with periodic Gaussian and 
approximating the Wilson model with arbitrary precision, does 
not exist. 



TABLE IV: Results for the coefficients /o (slope of the scaling 
curve) and /i (the intercept of the scaling curve) in the weak- 
coupling limit for the antisymmetric mass gap. 



Source 


fo 


A 


Villain (Hamiltonian) [3^ 


6.345 


4.369 


Morningstar [5Lj 


5.23 


5.94 


Hamer and Irving [59J 


5.30 


6.15 


Hamer, Oitmaa and Weihong 160] 


5.42 


6.27 


Plaquette Expansion ^lj 


5.01 


5.82 


Heys and Stump ]62J 


4.97 


6.21 


Lana [63] 


4.10 


4.98 


Xiyan, Jinming and Shuohong [64] 


5.0 


5.90 


Dabringhaus, Ristig and Clark 165] 


4.80 


6.26 


Darooneh and Modarres [66] 


4.40 


5.78 


Present Work 


5.39(9) 


6.3(1) 



practice, this limiting value is found by increasing [3 from 
strong-coupling until the mass ratio levels off in the weak 
coupling region. Earlier studies of the photon mass [fji| 
showed that the scaling of the mass ratio sets in for (3 > 1. 
If the continuum theory admits a stable bound state of 
two photons (a glueball), then the weak-coupling limit 
of the mass ratio will lie between 1 and 2. If the con- 
tinuum theory is simply a free-field theory of massive 
scalar photons, as in the Villain model or if the glueball 
remains in the continuum theory only as a resonance, 
then the ratio Rm — 2 should be observed as (3 becomes 
large. Lana |H3] argued for a sharp "transition point" at 
(3 = 1.40 where the symmetric bound state ceases to be 
stable, and crosses the level consisting of two free axial 
particles, so that mass ratio suddenly levels out at value 
2. 

In general, the most accurately calculated physical 
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quantity is the string tension K. Therefore the first quan- 
tity that we shall extrapolate to the continuum limit will 
be of the form M/K, where M is the glueball mass. The 
leading corrections to such a ratio are known to be of the 
order 0(1/ Q 2 j49|, where £ is some length scale. In the 
weak-coupling limit, 

M A /K^ir 2 /2 as (3 — > 00, 

in constrast to QCD, where M/y/K is constant in the 
continuum limit. Fig. [21 displays the estimates for this 
quantity, as a function of /?. The data certainly appear 
to level out below a 2 ,** ~ 0.006 at a value 6. This is quite 
close to the predicted value for the Villain model. 
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FIG. 14: The dimensionless ratio Ma/ (3K as a function of ef- 
fective lattice spacing, a e //- Open triangles show our present 
PIMC estimates. Our previous estimates using the Wilson 
action [l3| are shown as solid triangles. 



TABLE V: Comparison of the Hamiltonian estimates for di- 
mensionless mass ratio Rm obtained using the improved and 
standard Wislon actions 






Rm 

Improved action 


Wilson action 


1.0 


1.11(2) 


1.32(1) 


1.35 


1.4(1) 


1.5(2) 


1.55 


1.7(1) 


2.1(1.1) 


1.70 


1.9(2) 


2.2(1.5) 


1.90 


2.0(3) 


2.2(1.5) 


2.0 


2.0(5) 


2.4(1.6) 
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FIG. 15: Mass ratio, Rm, as a function of effective lattice 
spacing, a e //. Our present estimates are shown as the open 
triangles. The solid line is the fit to the data in the range 
0.00026 < a 2 e ff < 0.0066. The solid triangles show our earlier 
estimates obtained from the Wilson action [T3 |. 



Figure 1151 shows the behaviour of the dimensionless 
mass ratio Rm = Ms /Ma as a function of effective lattice 
spacing squared, CO/' wnere a eff is defined from Eq. 
as H3 

a eff = v/8^3e-" 2/3l,(0) . 

The plot shows that the mass ratio approaches very 
closely to the expected value of 2 in the large (3 limit. 
A linear fit to the data from 0.00026 < a 2 eff < 0.0066 
gives Rm = 2.007 ± 0.01, which is consistent with a con- 
tinuum limit value Rm = 2. There is no sign of a stable, 
scalar glueball bound state, nor of any sharp " break" in 
the mass ratio. In comparison with previous results using 
the unimproved action, there are two notable features. 
First, the results are much more accurate. Secondly, the 
corrections to the continuum limit appear to be linear in 
a%tf, whereas those for the unimproved action were lin- 
ear in a e ff This provides impressive evidence that 
the improved action has achieved what it was designed to 
do: produce faster convergence to the continuum limit. 



V. SUMMARY AND CONCLUSIONS 

In this work, we have applied the standard Euclidean 
path integral Monte Carlo method to obtain results in the 
Hamiltonian limit of the Symanzik improved U (1)2+1 lat- 
tice gauge theory on anisotropic lattices. Monte Carlo re- 
sults were obtained for the static quark potential, renor- 
malized anisotropy, the string tension, and the lowest- 
lying glueball masses. The inter-quark potential with 
the improved action exhibits good rotational symmetry. 
We found that both the improved and standard Wilson 
action do equally well in extracting the static potential 
at large separations. The improved discretization allows 
substantially more accurate estimates of the string ten- 
sion and the glueball masses. The extrapolations to the 
Hamiltonian limit were performed by simple quadratic 
fits. 

In this limit the string tension displays an asymptotic 
behaviour which is in excellent agreement with the be- 
haviour predicted by the weak-coupling approximation of 
Gopfert and Mack [33 • The scaling coefficient of the seal- 
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TABLE VI: PIMC results for the string tension, K, symmet- 
ric and antisymmetric glueball masses Ms, Ma and the mass 
ratio Rm- 



p 




K 


M s /K 


Ma/ K 


Rm 


1.0 


0.3724 


0.329(2) 


5.10(5) 


4.58(4) 


1.11(2) 


1.15 


0.2481 


0.189(3) 


7.8(1) 


6.1(1) 


1.28(3) 


1.25 


0.1884 


0.141(2) 


8.9(2) 


6.4(1) 


1.38(4) 


1.35 


0.1425 


0.112(5) 


10.8(5) 


7.3(4) 


1.4(1) 


1.40 


0.1239 


0.096(5) 


10.4(5) 


6.7(3) 


1.6(1) 


1.45 


0.1076 


0.067(4) 


13.9(9) 


8.0(5) 


1.7(2) 


1.55 


0.0810 


0.045(3) 


16.6(1.0) 


9.8(6) 


1.7(1) 


1.60 


0.0702 


0.041(2) 


18.1(1.1) 


10.1(6) 


1.8(1) 


1.65 


0.0608 


0.034(2) 


19.3(1.3) 


10.4(9) 


1.8(2) 


1.70 


0.0527 


0.028(2) 


24.1(2.1) 


12.6(1.0) 


1.9(2) 


1.75 


0.0425 


0.023(2) 


22.0(2.2) 


11.6(1.1) 


1.9(3) 


1.80 


0.0395 


0.021(2) 


21.6(2.7) 


11.1(1.3) 


1.9(3) 


1.85 


0.0341 


0.016(2) 


23.4(2.5) 


12.3(1.3) 


1.9(3) 


1.90 


0.0295 


0.013(2) 


25.9(3.2) 


13.0(1.6) 


1.6(3) 


2.0 


0.0220 


0.010(2) 


26.0(4.5) 


12.9(2.2) 


2.0(5) 


2.10 


0.0164 


0.009(2) 


23.4(5.0) 


11.5(2.5) 


2.0(6) 


2.20 


0.0122 


0.008(2) 


18.9(8.7) 


9.6(2.4) 





ing curve for the mass gap was estimated as roughly twice 
larger in magnitude than the weak-coupling prediction. 
This is an improvement over our previous unimproved re- 
sults which were estimated larger by a factor of 5-6. We 
believe that these estimates can be further improved by 
taking the renormalization of the couplings into account, 
i.e, computing the action beyond the tree-level. 

The weak-coupling behaviour of the antisymmetric 
mass gap was found to agree with the theoretical expec- 
tations. The parameters agree quite well with the results 
obtained from previous numerical and analytic calcula- 



tions. The mass ratio Rm = Ms /Ma was observed to 
scale to 2.007 ± 0.01 in the continuum limit. No sign 
of a glueball bound state or a sharp crossover between 
the levels was seen. This is in excellent agreement with 
the statement of Gopfert and Mack j3?J , that the contin- 
uum limit corresponds to a theory of free bosons, where 
Rm = 2 exactly. It also shows very clearly the value of 
using an improved action, in giving more rapid conver- 
gence to the continuum limit, as well as improved accu- 
racy. 

Taking advantage of the improved discretization on 
anisotropic lattices, we aimed to apply the Euclidean 
path integral Monte Carlo approach to examine the 
Hamiltonian limit of U(l) theory. The results obtained 
here clearly demonstrate that PIMC is a more reliable 
technique than other quantum Monte Carlo methods, 
such as Green's Function Monte Carlo, which gave only 
qualitative estimates of the string tension and the mass 
gap. 

Although suffering from the disadvantage that it re- 
quires an extrapolation to the anisotropic limit, PIMC 
is the preferred Monte Carlo technique for obtaining re- 
liable results in the Hamiltonian formulation, just as in 
the Euclidean formulation. In order to make the PIMC 
method a more valuable tool in Hamiltonian lattice gauge 
theories, it will be crucial to show that it allows one to 
treat eventually matter fields as well as gauge fields es- 
pecially in the non-Abelian case. 
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